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Abstract 

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 /im bands taken with the 
Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining 
the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on 
the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse 
depths of 0.197% ± 0.028% and 0.227% ± 0.025% in the 3.6 /xm and 4.5 fim bands, respectively. Our 
observations are best matched by models showing a hot dayside and, depending on our choice of 
model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted 
center of eclipse of 0.078 ± 0.032 hours, translating to e coscj = 0.0031 ± 0.0013 and consistent with a 
circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a 
previous transit study. 

Subject headings: eclipses — planetary systems — WASP-5b — techniques: photometric 



1. INTRODUCTION 

Hot Jupiters are a class of extrasolar planet that, 
as the name suggests, are similar in size and composi- 
tion to Jupiter but orbit very close to their parent star 
and have correspondingly high effective temperatures, 
ranging from 1000K to over 3000K. By measuring the 
wavelength-dependent decrease in light during the sec- 
ondary eclipse (when the planet passes behind its star), 
we can characterize the planet's em ission spectrum and 
deduce its atmospheric p roperties (jDeming et al.l [20051 : 
ICharbonneau et al.ll2005l ). These atmospheric properties 
include: whether or not the atmosphere has a tempera- 
ture inversion and how well heat is redistributed from the 
planet's dayside to its nightside. In addition, the relative 
timing of transits and secondary eclipses can constrain 
the eccentricity of the planet's orbit. 

During its cryogenic mission, Spitzer observed sec- 
ondary eclipses for 15 extrasolar planets. The re- 
sulting studies indicate that hot Jupiter atmospheres 
can be characterized by the presence or absence of 
a strong thermal inversion in the upper atmosphere 
(IBurrows et al.ll200aiFortnev et al.| [2008t IBarman 112001 
Madh usudhan et al.l I2009TT Although it exhausted the 
last of its cryogen in 2009, the Spitzer Space Telescope's 
Infrared Array Camera (IRAC) remains f unctional in its 
3.6 and 4.5 /xm bands ([Fazio et al.l [2004) . and the tele- 



scope has continued to survey the emission spectra of hot 
Jupiters as part of Spitzer's extended warm mission. 

Some groups have suggested that absorbers such as 
TiO in the upper atmosphere of hot Jupiters are re- 
sponsible for atmosp heric inversions ([Hubenv et al.|[200l 
iFortnev et al.l [20081 ). However, large-scale atmospheric 
mixing would be required to preserve gaseous TiO in 
the upper atmosphere, as this molecule should condense 
in the deep interiors of most hot Jupiters. It is uncer- 
tain whether such macroscopic mi xing should take place 
in a stably stratified atmosphere dShowman et all 120091 : 
ISpiegel et al. l[2009l : IParmentier et al.ll2013D . In addit ion , 
temperatur e inversions have bee n observed on planets 
like XO-lb (jMachalek et al.ll2008f ) that have dayside tem- 
peratures below the condensation point of TiO. Sulfur- 
containing comp ounds have been pr oposed as an alter- 
native absorber (jZahnle et al. Il2009| ). 

In this paper, we present observations of the transit- 
ing hot Jupiter WASP-5b. This planet is very dense 
compared to other planets in its class, suggesting the 
presence of a large, metal-rich core (|Anderson et alj 
2008t iSouthworth et al.l 120091: iMadhusudhan et al.ll2009l: 
Gillon et all 120091: iDragomir et al.l 120111 iF ukui et al l 
20111: iPont et all 120111 : iHover et all 120121 ) . Knutson et 
al. 2010 reported evidence for a correlation between 
stellar activity and hot Jupiter emission spectra, where 
hot Jupiters with strong temperature inversions tend to 



2 



TABLE 1 

Best-fit values for the secondary eclipse depth and time. 



Wavelength (^tm) 


Eclipse Depth (%) 


Brightness Temperature (K) 


Center of Eclipse (BJD) 


Eclipse Offset (hours) 


3.6 


0.197 ±0.028 


2210 ± 140 


2455200.42472 ± 0.0021 


0.077 ±0.039 


4.5 


0.227 ±0.025 


2090 ± 120 


2455174.36986 ± 0.0021 


0.080 ± 0.060 
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Fig. 1. — Raw photometry at 3.6 and 4.5 fim vs time from 
the center of the predicted secondary eclipse. The decorrelation 
functions to correct for intrapixel sensitivity are overplotted in red. 
Data are binned in sets of 25 points. 

orbit more quiet stars and more active stars typically 
host planets without inversions. This may be due to 
the fact that the more intense UV radiation from ac- 
tive stars destroys the compounds that are responsible 
for creating thermal inversions. WASP-5 is a modestly 
active G4V star with log(R') = -4.75 and Shk = 0.215, 
so it is not expected that this planet would have a 
strong tempera t ure in v ersion. In addition, papers s uch as 
Cowan fc Ago! (l20ll . IShowman fc Guillot I l|2002T) . and 
Perna et al.l (|2012f ) find that hotter planets have weak 
transport of energy to the nightside. With an effective 
temperature of 1717 K, we would expect WASP-5b's at- 
mospheric energy transport to be lower. In this paper, 
we wi ll tes t the correlations prop osed by iKnutson et al.l 
(2010) and lCowan fc Ago] (|2011h by constraining the at- 
mospheric properties of WASP-5b . 

In a previous study, iFukui et"aTI (j2011l ) found that the 
intervals between WASP-5b's transits do not appear to 
be constant, which could indicate the presence of an ad- 
ditional body in the system perturbing the planet's or- 
bit. Hoyer et al. (2012), however, disputes the claims 
of such transit timing variations (TTVs) . In Section 13.21 
we present the first measurements of this planet's sec- 
ondary eclipse times and discuss the implications of our 
measurements for the proposed perturber. 

2. OBSERVATIONS AND METHODS 

We began by extracting photometry from calibrated 
images in the 3.6 fim and 4.5 /im bands from the Spitzer 
Space Telescope. Each channel consists of 2115 images 
acquired over a period of 7.7 hr. The basic calibrated 
data (BCD) files used were dark-subtracted, linearized, 
flat-fielded, and flux-calibrated. We correct for transient 
"hot pixels" in a 20 x 20 pixel box around the star by 
removing intensity values > 3er from the median value of 
the surrounding frames. These values are then replaced 
by the median. In addition, to further reduce the noise in 
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Fig. 2. — Photometry and best-fit eclipse curve for both wave- 
bands after decorrelation vs time from the predicted center of the 
secondary eclipse. The data are binned in sets of 25 points. 

the data we trimmed the first 53 minutes from channel 1 
data and the first 15 minutes from channel 2 data, which 
exhibit larger deviations in position. 

We performed aperture photometry with radii rang- 
ing from 2.0 to 3.0 pixels in intervals of 0.1 pixels, and 
from 3.0 to 5.0 pixels in half-pixel intervals. We found 
that setting the aperture for channel 2 equal to 2.2 pixels 
produced the smallest root mean square (rms) scatter in 
the data. We found that we obtained superior results in 
channel 1 using a time- varying aperture proportional to 
the square root of the n oise pixel value (jKnutson et al.l 
[20H iLewTs et al. \\20l3 \ minus 0.1 pixels. In channel 
2, however, using a time varying aperture increased the 
scatter, so a fixed size aperture was selected. Our me- 
dian aperture size for channel 1 was 1.7 pixels. We de- 
termined the position of the star in both channels using 
flux- weighted centroiding with a radius of 5.0 pixels. 

We next correct for the dominant instrument effect in 
the two bands. Since the sensitivity of individual pix- 
els varies from the center out to the edge, the apparent 
flux of the star will fluctuate with the movement of the 
telescope. To correct for this intrapixel sensitivity varia- 
tion, we fit the data with quadratic functions in x and y 
position. For channel 1, the measured flux / is given by 

/ = fo (ci (x ~ x ) 2 + c 2 (y - yo) 2 + c a (x - x ) + c 4 (y -yo) + cs) 

where fo is the incident flux, x and y are the positions of 
the star on the array, xq and yo are the median positions 
over the time series, and ci_5 are free parameters. For 
channel 2, we found that adding a factor (x — xo){y — yo) 
further reduced the rms scatter in the data. We also 
tried fitting cubic and linear functions, but a quadratic 
function yielded the best fit. Figure [1] shows the raw data 
together with our best-fit quadratic functions. Figure [5] 
displays the normalized data with best-fit eclipse models 
for each channel. 

We us ed a Mark o v Chain Monte Carlo (MCMC) 
method (|Ford1 [20051 iWinn et all l200l with 10 5 steps 
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Fig. 3. — Dayside planet/star flux ratio vs wavelength for four 
model atmospheres (Fortney et al. 2008) with the band-averaged 
flux ratios for each model superposed (colored circles). The ob- 
served contrast ratios are overplotted as the black circles, with 
uncertainties shown. Two models contain gas- phase TIP, which 
produ ces a high-altitude temperature inversion (Parmentier et al. 
120131) . In the figure, the red and blue models represent a planet 
that radiates heat only from the dayside. The red and blue models 
represent a planet that radiates heat only from the dayside, while 
the green and purple models distribute heat equally between the 
two hemispheres. 

to determine the best-fit parameters and uncertainties. 
We used seven free parameters for channel 1, and eight 
parameters for channel 2, including the eclipse depths, 
timing offsets, and the intrapixel sensitivity corrections. 
The planetary and stellar radii, orbital period, and or- 
bital inclination were set to the values given in Hoyer 
et al. (2012). We took the median value yielded by 
the MCMC as the best-fit solution, and the symmetric 
range over which the probability distribution contained 
68% of the points around the median as the uncertainty. 
The distributions for all parameters were approximately 
Gaussian. See Table Q] for the best- fit eclipse depths and 
times. 

We also calculate error bars using the "prayer-bead" 
method (Gillon et al. 2009), which involves shifting the 
residuals of the data in single point increments and recal- 
culating the best-fit eclipse depth and time. This method 
yields error bars consistent with those of the MCMC, and 
we chose to use the larger of the two errors in each case. 
For channel 1, we used the prayer-bead error for both 
the eclipse depth and offset (0.028% and 0.039 hours, re- 
spectively) as opposed to the MCMC error (0.027% and 
0.038 hours). For channel 2, we used the MCMC for the 
eclipse depth (0.025% as opposed to 0.024% from the 
prayer-bead) and for the eclipse offset (0.060 hours as 
opposed to 0.045 hours from the prayer-bead method). 
The close agreement between these two methods indi- 
cates that there is minimal time-correlated noise in our 
data. 

3. DISCUSSION 

3.1. Atmospheric Temperature Structure 

We consider two classes of atmospheric models in this 
paper. The first, following Fortney et al. (2008) exam- 
ines the impact of TiO as an absorber in the upper atmo- 
sphere. These models parameterize the redistribution of 
heat by adjusting the incident flux at the top of the at- 
mosphere. In Figure |3l four models are shown. The red 
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Fig. 4. — Dayside planet/star flux ratio vs wavelength for three 
model atmospheres (Burrows et al. 2008) with the band-averaged 
flux ratios for each model superposed (colored circles). The ob- 
served contrast ratios are overplotted as the black circles, with 
uncertainties shown. The model parameter k is related to the at- 
mosphere's opacity, while p is related to the heat redistribution 
between the day and night sides of the planet (p = 0.0 indicates no 
heat redistribution, and p = 0.5 indicates complete redistribution) 

and blue models represent a planet with a hot dayside 
that redistributes little heat to the nightside. The green 
and purple models characterize an atmosphere in which 
heat is evenly distributed between the planet's day and 
night sides. The blue and purple models contain gas- 
phase TiO, which acts as an absorber in the planet's 
upper atmosphere. The red and green models contain no 
TiO. We find that the slope between the band averaged 
flux ratios on the blue model is too steep to fit both of 
our observed values. As seen in the figure, our observa- 
tions are best matched by the red model, which has no 
thermal inversion and poor redistribution of energy from 
the dayside to the nightside. 

The second class of models follow the methods of 
iBurrows et al.l (|2008[ ). As opposed to modifying the 
atmospheric levels of TiO, this method utilizes a non- 
specific gray absorber at low pressures. A heat sink is 
added at depth to allow for the redistribution of heat 
from the planet's dayside to its nightside. We use the di- 
mensionless parameter P n to specify the heat redistribu- 
tion, where P n = 0.5 indicates evenly distributed energy 
among both the day and night sides of the planet, and 
P n = 0.0 indicates no heat redistribution. This model 
also puts an extra absorber with an optical opacity, n, 
in the upper atmosphere to heat it by partial absorption 
of the incident stellar fiuxin visible wavelengths. An en- 
hancement in optical wavelength opacity in the planet's 
atmosphere will produce a thermal inversion, but impor- 
tantly raises the temperature of the outer atmosphere 
to better reproduce the photometric signatures of a sub- 
set of measured hot Jupiters. The character of such an 
absorber, and whether it is a molecule or haze, are cur- 
rently unknown.. Figure 3] shows three of these models. 
Our observations are best matched by the green model, 
which indicates a weak temperature inversion and little 
heat distribution be tween the sides of the planet, in good 
agreement with the lFortnev et all (|2008l ) models. 

We calculate brightness temperatures for the planet 
in both bands using a PHOENIX mode l atmosphere 
for the star with stellar parameters from IDovle et al. I 
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(2012) and find planetary brightness temperatures of 
2210 ± 140 K at 3.6 fim and 2090 ± 120 K at 4.5 /im. 
If the planet has an albedo of zero and re-radiates the 
absorbed flux uniformly over its entire surface, we cal- 
culate a predicted equilibrium temperature of 1700 K. 
If it re-radiates this energy from the dayside hemisphere 
alone, this temperature increases to 2030 K. The mea- 
sured brightness temperatures therefore provide a model- 
independent indication that WASP-5b has a hot dayside 
and relatively inefficient day-night recirculation. We con- 
sider whether or not this is typical of the class of hot 
Jupiters by plotting the measured surface fluxes for all 
planets with published secondary eclipses in the 3.6 and 
4.5 fim Spitzer bands and comparing these fluxes to the 
values for a blackbody at the predicted equilibrium tem- 
perature (see Figure [5]). We find that these planets uni- 
versally exhibit flux excesses relative to the coolest pre- 
dicted blackbody (corresponding to efficient day-night 
circulation) in at least one, and sometimes both, Spitzer 
bands. We do not find any clear correlation between flux 
excesses and the predicted planet temperature, although 
planets orbiting active stars tend to lie towards the lower 
right region of the plot. 

Our results are consistent with the correlation between 
stellar ac tivity and hot Jupiter emission spectra, as de- 
scribed in lKnutson et al.l (|2010|) . According to this study, 
hot Jupiters with strong temperature inversions tend to 
orbit more quiet stars, while more active stars typically 
host planets without inversions. WASP-5 is a G4V star 
with log(R') = —4.75 and Shk — 0.215, indicating that 
the star is modestly active. This is consistent with our 
results from the atmospheric models (Section l3.1[) . which 
suggest the presence of either a weak thermal inversion 
or none at all. In addition, the low heat redistribu- 
tion of WASP-5b is consistent with the predictions of 
ICowan fc Ago! (pOll . 

3.2. Orbital Eccentricity 

The timing of the secondary eclipse depends on the ec- 
centricity of the planet's orbit. When accounting for the 
additional 26.78 s that l ight takes to travel across WASP- 
5b's orbit (|Loeb Il2005l ), a secondary eclipse for circular 
orbit would occur at a phase of 0.5002. As a result, we 
can utilize the predicted and observed times of eclipse 
to constrain ecosoj, where e is the orbital eccentricity 
and oj is the argument of pericenter (Charbonneau et al. 
2005). 

We observed a mean offset from the predicted cen- 
ter of eclipse of 0.078 ± 0.032 hours, translating to an 
ecosw value of 0.0031 ± 0.0013. Despite ongoing circu- 
larization from tidal interactions with the parent star, 
the gravitational influence of a second planet could pre- 
serve a non-zero eccentricity. The presence of such tran- 
sit timing variations due to an additional perturber in 
the system has been disputed in past articles (Fukui et 
al. 2011, Hove r et a l. 2012). Following t he model of 
ICowan fc Agoll (|2011l ). as implemented in lAgol et al. I 
(|2010D . we calculate that the maximum expected time 
offset due to an offset hot spot would be 37 s or less than 
1 sigma for our measurement. We conclude that ecosw 
is consistent with zero. This indicates that, unless our 
line of sight is aligned with the planet's semi-major axis, 
the orbit is circular. 
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Fig. 5. — Ratio of measured to predicted planet surface flux in 
the 3.6 vs. 4.5 fimSpitzer bands for all planets with published sec- 
ondary eclipse detections at these two wavelengths. The colors of 
the points reflect the predicted equilibrium planet temperature T eq , 
which is calculated assuming that the planet has an albedo of zero 
and re-radiates energy uniformly over its entire surface. Predicted 
fluxes are calculated assuming the planet radiates as a blackbody 
at this temperature. Measured fluxes are calculated using pub- 
lished secondary eclipse depths and assuming a PHOENIX atmo- 
sphere model for the star. For planets with eccentricities greater 
than 0.05 measured at 3cr significance or higher we show the orbit- 
averaged temperat ure. Planets orbitin g active stars (defined as 
having \og(R'^^) (Knutson et al. 2010) greater than —4.9 and a 
stellar effective temperature less than 6200 K) are shown as filled 
squares, all other planets are plotted as filled circles. 

4. CONCLUSIONS 

We observed secondary eclipses of the extrasolar planet 
WASP-5b in the 3.6 and 4.5 /an bands, and compared 
our measurements with atmospheric models generated 
according to the methodology of Burrows et al. 2008 
and Fortney et al. 2008. Our eclipse values of 0.197% ± 
0.028% and 0.227% ± 0.025% in these bands are best 
matched by models indicating a weak thermal inver- 
sion or no inversion at all. Our measurements are 
consistent with the observed correlation between stel- 
lar activi ty and the presence o f thermal inversions pre- 
sented in iKnutson et al.l (|2010| ). In addition, these one- 
dimensional atmospheric models all indicate a hotter 
dayside temperature for WASP-5b, with the Burrows 
models suggesting a redist ribution factor o f 0.1, consis- 
tent with the prediction of ICowan fc Ago] (|2011f ). Mea- 
surements at additional wavelengths would be helpful to 
confirm or reject the presence of the inversion. These re- 
sults further demonstrate the utility of warm Spitzer in 
characterizing the properties of hot Jupiter atmospheres, 
despite being limited to the 3.6 and 4.5 channels. Ad- 
ditional measurements in the J, H, and K bands would 
also allow us to constrain the metallicity of WASP-5b's 
atmosphere, which is likely to be enhanced relative to 
other hot Jupiters. 

By measuring the mean timing offset from the pre- 
dicted center of eclipse (0.078 ± 0.032 hours), we calcu- 
late the parameter ecosw to be 0.0031 ± 0.0013. This 
indicates that the planet's orbit must be nearly circu- 
lar unless our line of sight is well-aligned with the major 
axis,. These observations, along with those in Hoyer et 
al. 2012, argue against the possibility of ongoing pertur- 
bations from an additional body in the WASP-5 system. 
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